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A geometrical optics model is presented that is based on defining a 
local numerical aperture as a junction of fiber radius with a uniform 
power distribution. Transmission vs transverse offset characteristics for 
parabolic-profile fiber splices are calculated for unequal fiber-core diame- 
ters. We show that the often-used assumptions of equal-mode excitation, 
equal-mode attenuation, and no-mode coupling are not adequate to calcu- 
late realistic transmission vs offset characteristics. Splice-loss measure- 
ments with long fiber lengths on each side of the splice show less than the 
calculated sensitivity to small offsets and greater than the calculated sensi- 
tivity to large offsets. 

I. INTRODUCTION 

Accurate transverse alignment is difficult to obtain in fiber optic 
splices primarily due to the small size of optical fibers. Many splicing 
techniques 1 - 2,3 use grooves or channels and are dependent on fiber outer 
diameter (od) for transverse alignment. Experience thus far in- 
dicates that od variations as small as ±1 percent can be achieved on 
a single long fiber; however, fiber-to-fiber variations for fibers drawn 
at different times may be considerably greater. Small fluctuations in 
core-diameter-to-OD ratios are also expected. When the receiving 
fiber core is smaller than the transmitting fiber core, a loss occurs even 
with perfect transverse axial alignment. 

Thiel 4 and Henderson 5 have reported transmission vs offset for equal- 
core-diameter, step-profile fiber splices and Pugh 6 has considered core- 
diameter ratios less than one (receiving fiber smaller than transmitting 
fiber) for step-profile fiber splices. Calculated values for transmission 
vs offset for equal-core-diameter, parabolic-profile fiber splices have 
been obtained earlier by the author. 7 
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Transmission vs offset for parabolic-profile fiber splices with unequal 
core diameters is calculated in this paper. In addition, transmission vs 
offset-measurement data for equal-core-diameter, parabolic-profile 
fiber splices are presented. We found that the theoretical model gives 
pessimistic results for the equal-diameter case for the region of primary 
interest (offsets less than 0.8 core radius). Measured transmission is 
less than the model for offsets greater than 0.8 core radius. 

II. ASSUMPTIONS 

The solid angle defined by the local numerical aperture (na) of the 
fiber at every point on the core is assumed to contain a uniform power 
distribution. This assumption is consistent with equal-mode excitation, 
equal-mode attenuation, and no-mode coupling. 8 Consider the profile 



(1) 



where 



n(r) = mU ~ 2A (^) a ] for r<R ' 

A = (vo — Ve)/vo is small. 
ilo = refractive index at center of core. 
77 c = refractive index of cladding. 
a = a parameter between 1 and °o . 
R = fiber core radius. 

The na as a function of radius is then 

NA(r) - [,<r)« - ift]*, (2) 

or 

xt(r) £^oV2A[l- (g)"]*' (3) 

Since uniform power per unit solid angle is assumed transmitted for 
each incremental area of the core within the angle 

0( r ) - sin-* [NA(r)A(f)], (4) 

the total power, Pt, equals 

Prsp(o) /,7'[ i -(s)1'**' (5) 

or 

P T = xfl 2 aP(0)/(a + 2), (6) 

where P(0) is a constant dependent on the input power A and ij . 
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When the fiber cores arc offset, the total power received at a given 
point is assumed to be limited by the transmitting or receiving na, 
whichever is the minimum at that point. 

III. GENERAL APPROACH FOR THE PARABOLIC PROFILE (a = 2) 

For the case of equal fiber-core diameters, the locus of equal na is a 
straight line, as shown in Fig. 1. The na function is integrated over the 
area of overlap bounded by a circle and a straight line. The resulting 
integral for received power in region I, where the transmitting fiber has 
the minimum na, is 



/•cos- 1 (dttRr) rRT r / r \ 21 

JO J d/2cosfl L V Rt / J 



rdrdd. 



(7) 



For region II, where the receiving fiber na function has the minimum 
na, a positive translation of the na function must be made. The re- 
ceiving na function, translated an amount d and referenced to the 
center of the transmitting fiber, equals 



NA(r) ^n V2A |l - ( 



r 2 - 2dr cos B + d 2 



)i 



(8) 



Four separate cases are considered to calculate transmission vs off- 
set for unequal core diameters. If K is the ratio of the radius of the 
receiving fiber to the radius of the transmitting fiber, K = R R /R T , 
and d = offset, then the four cases are as follows : 



TRANSMITTING 
FIBER CORE 




RECEIVING 
FIBER CORE 



Fig. 1 — Regions of overlap for offset fiber cores (Rt = Rr). 
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Case I 

See Fig. 2. d/R T < K — 1 f or K > 1. Core boundaries do not inter- 
sect and Rr > Rt- 

Case II 

See Fig. 3. d/R T > K — 1 for K > 1. Core boundaries intersect 
and Rr > Rt. 

Case III 

See Fig. 4. d/Rr < I — K for K < 1. Core boundaries do not inter- 
sect and Rt > Rr- 

Case IV 

See Fig. 5. d/Rr > 1 — K for K < 1. Core boundaries intersect 
and R T > Rr. 

For Cases I and III, areas of overlap are circles and the integrals 
are easily written and solved for the transmission. The resulting equa- 
tions are shown on Figs. 2 and 4. For Cases II and IV, some areas of 
overlap are bounded by arcs of circles, and integrals for these areas 
are tedious to solve. The approach used for these cases is to divide 



K= R R /R T 
a = d/(1-K) 
w = d/(1 +K) 
u = (a - w)/2 
q = w- u 



LOCUS OF- 
EQUAL NA 



*-#-(# 




,2 2<d-u)rcos£ 





2(^+_P r 

P TOT 



R R H R T R T 



T - fO'. i + is. 17, _ 1 )(i + 5au + "fu ± fia. + u _ i\\ 

Fig. 2— d < Rr - R T . 
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k= R R /R T 

u = (1/2d)(R T 2 -R R 2 +d 2 ) 

w = d/(1 + K) 



2uw - w 2 - R 2 

C= 2(u-w) 
p = c — w + u 



Fig. 3— d > Rr- R t . 



the areas of overlap into sections which are bounded by a circle and a 
straight line and to combine these sections to obtain the required 
geometry. All integrals then have the form suggested by eqs. (7) and 
(8) and can be segmented into the following three integral forms. 



a a 



rdrdd = ir cos" 1 j--^ V6 2 - a 2 . 

i/cosS ^ 1 



/•coa * (a/6) r t 

J a/ 

/■cos -1 (a/ 6) f b 

/ / r 2 cos 6 drdQ = \ V(6 2 - a 2 ) 3 . 

rcoi-Ualb, rb & 4 „ a 

/ / r 3 drdd = -r- cos -1 r - ts V (6 2 — a 2 ) 3 

JO J a/cos e LZ 



4 



(9) 
(10) 

(11) 



These integrals were computer-programmed and summed to calcu- 
late transmission vs offset. Closed form expressions were not obtained 
due to the large number of terms involved. The combination of in- 
tegrals for cases II and IV is given in the Appendix. 
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V*//[l-(*^^-£H« •• P n = 2//[l-(^ R ) 2 ]rdrd fl - 2 //[,-(y : 



+ 2(»-d>rco,g _(u^d}'i fdrdfl 
Rr 2 r r J 

Fig. 4— d < i2 r - fl«. 

IV. RESULTS OF CALCULATIONS 

Figure 6 is the resulting plot of transmission (percent) vs offset 
(transmitting fiber-core radius) for K = 1.05, 1.00, and 0.95. Also 
shown are results calculated for the step profile for K = 1.00 and 0.95. 
These calculations show a 33-percent greater sensitivity to offset for 
the parabolic profile and approximately the same sensitivity to K as 
the step profile. Figure 7 is a family of curves for small offsets. For zero 
offset, the transmission equals 1 for K > 1 and equals K 2 for K < 1, 
as in the step profile. For offsets less than Rt — Rr for K < 1, the 
parabolic profile exhibits a dependence on offset, whereas the step 
profile is constant. The reduction in transmission due to this effect is 
only 0.7 percent at 0.05 core-radius offset for K = 0.95. The slope 
discontinuities in the parabolic case in Fig. 7 result from the 0.01 core- 
radius offset increment selected. Curves for transmission vs offset are 
smooth except at the point d = Rt — Rr, where a discontinuity in the 
slope may exist. 
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k = R R /R T 

u = (1/2d)(RT -Rp +d 2 ) 

w = d/(1 + K) 



-2uw + w 2 + R^ 
C 2(u-w) 

p = c + w — u 



Fig. 5— d >R T - Rr. 



PARABOLIC 



O O O STEP 




Fig. 6 — Calculated transmission vs transverse offset. 
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Fig. 7 — Calculated transmission vs transverse offset for small offsets. 

V. COMPARISON OF CALCULATIONS AND MEASUREMENTS 

Several investigators have seen indications 9 ' 10 that the sensitivity 
to small offsets is significantly less than the calculated values for para- 
bolic fiber splices. These measurements 9 - 10 were made with short fibers 
on each side of the splice and the launched mode distribution probably 
differed greatly from the uniform power distribution assumed here. 

A measurement of transmission vs offset was made with 500 meters 
of Corning fiber on each side of the splice subsequent to measuring the 
100-percent transmission level. A helium-neon laser with an expanded 
beam and a 20X objective lens were used to completely fill the fiber 
na. Two tabs were placed near the center of the 1-km parabolic cow 
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fiber, about 4 inches apart, and the fiber was scored and broken be- 
tween these tabs. A single break was made and the same ends were 
realigned on a microscope slide against a thin straightedge with the 
tabs used to obtain rotational alignment. Glycerin was used as the 
index-matching material and a cover slip held the ends in position 
under a microscope. The 100-percent level was again obtained and 
transverse offset was introduced in the splice and photographed. 

Figure 8 shows the result of the transmission vs offset measurement. 
The sensitivity to offset for small offsets was less than predicted by 
theory, and greater than predicted for large offsets. The crossover 
point is approximately 0.8 core radius. For small offsets the reduction 
in sensitivity to offset is approximately the same as one would expect 
for a 10-percent oversized receiving fiber (K = 1.1). 

The use of the geometrical optics approximation has been corrobo- 
rated by comparison with the wave analysis of Marcuse. 12 The con- 
clusion is that the assumption of equal-mode excitation, equal-mode 
attenuation, and no-mode coupling is not adequate to calculate splice 
transmission sensitivity to offset. The calculated transmission with 
core-diameter mismatch shows trends that can be used to obtain rela- 
tive effects only. The long fiber length measurements of transmission 
vs offset can be used as design data for setting required alignment 
tolerances. 
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Fig. 8 — Comparison of calculated and measured transmission vs transverse offset. 
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Current work by Gloge 13 indicates that the assumption of a sharp 
cutoff of power at the boundary of the solid angle defined by local na 
could cause the discrepancies noted here. The model used in this 
paper is extended by Gloge to include a nonuniform power distribution 
across the solid angle defined by the local na. 
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APPENDIX 

Integrals for Cases II and IV 

A.1 Case II (Fig. 3) d/Rr >1-Klor K>1. 

/•cos- 1 (-u/flr) rRT T / r \2~| 

/ h - ( p- ) rdrde 

/•co fl -i(p/e) re r ^ r 2 _ 2r(c — w) cos 6 

JO Jp/cobB L Rt 



R% 

(c- 



R% J 
roof 1 (ph) re r r 2 2r(c + d — w) cos 

Pll = 2 I / 1 — pa 

JO J p/coa 6 l_ B>R 



rdrdd 



«R 

(c + d- w) 2 ' 



rdrdd 



/•cob" 1 [ (d-u) I Rr) i-Rr r / r \ 2 T 

/ h - ( w ) rdrde 
D J (d-u) /cos 9 L \ *** / J 

T _ Pi + Pii + Pin 
Pt 

A.2 Case IV (Fig. 5) d/R T >1 -KforK<1. 

rcoB- 1 [(,u-d)IRR] rRRV / r \2"| 

- 2 /. /. t 1 - U) r r<M 

rcoB-Upie) re r / r \ 2 2r(c + w - d) 
io Jp/coae |_ \RrJ Rr 

(c + w- d) 2 
R% 

926 THE BELL SYSTEM TECHNICAL JOURNAL, SEPTEMBER 1976 



COS 



rdrdd 



/•cos l (pic) re r / r \ 2 2r(c + w) cos 

(c + wy 1 
B| J 

/■cos -1 (u/fir) rflr r / r \ ! 

Pl " - 2 / / x " ( RZ ) 



rdrdd 



rdrdd 



T _ Pi + Pii + Phi 
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